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Technical Work 

During this period efforts have been expended for the 
following tasks: 

(1) Microstrip Antenna Module 

In February, an initial computer program to be used in 
analyzing the four-element array module was completed. This 
program performs the analysis of modules composed of four 
rectangular patches which are corporately fed by a microstrip 
line network terminated in four identical load impedances. These 
identical load impedances are determined by ignoring mutual 
couplings which are believed to be second-order effects and by 
applying the cavity model theory to an individual, isolated 
patch. The resulting impedance at the module input is then 
determined by using transmission line theory and lumped, lossless 
circuits to represent the microstrip line corners and junctions. 

Currently, a rigorous full-wave analysis of various types of 
microstrip line feed structures and patches is being carried out. 
As shown in Figure 1.0, these include the microstrip line feed 
and the electromagnetically coupled microstrip line feed between 
layers of different electrical parameters. The configuration for 
the array module remains unchanged except that the method of 
moments is now used to determine the input impedance. Mutual 
coupling is ignored at the moment, although it could be taken 
into account with more computational effort, and TEM transmission 
line theory is still used. 
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In the method of moments analysis, the distribution and 
magnitude of a primary electric or magnetic current source are 
assumed as the antenna excitation (e.g., electric current 
filament for coaxial probe feed or magnetic current ribbon for 
narrow feed gap), and the equivalence theorem is used to replace 
the antenna and feed lines by an unknown secondary current 
distribution. This secondary current distribution, which 
accounts for the scattering of the fields of the primary source 
by the antenna and feed structure, will in general contain two 
orthogonal components of current which flow in the plane of the 
antenna. Consequently, a vector dyadic Green's function must be 
used to obtain the fields of these individual components. By 
imposing the required boundary condition that the total 
tangential electric field on the surface of the antenna and the 
feed line is zero, a set of integral equations relating the 
components of the fields of the primary and secondary sources is 
obtained. The solution of the integral equations then yields an 
approximation to the vector components of the total current 
flowing on the structure. Unfortunately though, the solution of 
these currents requires the evaluation of integrals containing 
integrands which are oscillatory and very slowly decaying. But 
once these currents have been obtained, the input impedance may 
then be determined since the source excitation was assumed to be 
known . 

At present, a method of moments program has been implemented 
for the case of a single dielectric layer and microstrip line fed 
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rectangular patches in which the primary source is assumed to be 
a magnetic current ribbon across the line some distance from the 
patch. It has been found that in order to obtain input impedance 
values which agree with measured values the location of this feed 
model on the feed line must approximately be at least one 
dielectric wavelength away from the patch edge. Initial 
experiments have been conducted to examine this problem and 
determine at what distance from the patch along the microstrip 
line the fields are approximately TEM. Typical results are shown 
in Figure 2.0, where it is seen that impedance values can be 
translated by simple TEM transmission line theory wth reasonably 
good accuracy. Further work on the theoretical modelling of the 
excitation is still required. 

A typical example comparison of measured values with those 
computed by the program is shown in Figure 3.0. In this program, 
use has been made of subsectional "roof-top" basis functions as 
shown in Figure 4.0; therefore, the program has also been used to 
examine a few types of microstrip discontinuities as well. A 
second program which analyzes multilayered dielectric structures 
of rectangular patches and electromagnetically coupled microstrip 
lines is nearly completed. Like the first completed program, 
this version also uses subsectional basis and testing functions 
and is completely self-contained, needing no external libraries 
to perform complex matrix inversions, numerical integrations, or 
computation of special functions. 



Both experimental and theoretical studies will continue 
concerning these feed structures. In particular, experimental 
work will be conducted to investigate theoretical primary source 
models other than the delta gap generator, the validity of the 
TEM transmission line model for different feed structures (such 
as co-planar waveguide) , and the validity of the assumption of 
insignificant mutual coupling in alternate array modules. 
Amendments to the programs themselves are intended, such as the 
use of entire domain basis functions to improve computational 


efficiency 
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Figure 1.0 Physical structures of different feed networks being analyzed. 

(a) Patch and feed networks are constructed on a single dielectric 
slab, (b) Patch and feed networks are constructed on two 
dielectric layers having different electrical parameters. 
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A- Measured at plane 
a- 3*2 cm, b* 2.15 
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Figure 2.0 Comparison of input impedance loci measured at planes A and B 
where the values at plane A have been translated to plane B by 
simple TEM transmission line theory, (a) Substrate is 0.034 X* 
thick at 4.0 GHz. (b) Substrate is 0.069 thick at 4.0 GHz. 
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Figure 4.0 (a) The basis function layout utilized in the method of moments 
analysis, (b) The tee junction and step in width types of 
microstrip line discontinuities which are analyzed through the 
use of subsectional basis functions. 
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(2) An Improved Theory for a Thick Circular Microstrip Antenna 
Fed with a Coaxial Cable 1 

CHAPTER 1 
INTRODUCTION 

The first twenty years, following the introduction of the microstrip antenna in the 
early 1950’s, were marked by few attempts to understand and analyze its characteristics. 
By the early 1970’s, however, technological progress in the fields of integrated microwave 
electronics and mm-wave engineering had created the need for a new generation of 
antenna systems. The microstrip antenna, because of such unique properties as light 
weight, low profile, ease of manufacture, compatibility with integrated microwave com- 
ponents, etc., was seen as a promising candidate to satisfy the emerging requirements. 
This realization gave impetus to the tremendous expansion of microstrip antenna research 
and development witnessed in the past fifteen years. During this time virtually every 
engineering aspect of this antenna problem has been explored to some degree. Neverthe- 
less, the major portion of the existing microstrip antenna literature is concerned with one 
specific topic, namely the input impedance problem. Interest in the input impedance prob- 
lem is stimulated both by the designer’s need for accurate characterization of the 
antenna’s terminal properties, as well as by the theoretical challenge that the problem 
presents to the antenna engineer. The reason for this challenge becomes clear when it is 
recalled that precise calculation of the input impedance requires the solution of a formid- 
able boundary value problem for the near fields of the antenna. Complicating the solu- 
tion even further is the fact that the value of the input impedance is critically dependent 
upon the somewhat arbitrary choice of the antenna’s input terminals. In fact [l] "to be 
able to define an antenna as a two-terminal element, it is necessary to define two electri- 
cally close points in the excitation region in which the feeder terminates and the antenna 
begins. This can be done fairly arbitrarily if rough analysis only is needed, but consider- 
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able care has to be exercised if more accurate analysis is required. 

In spite of the assumed electrically small volume occupied by the excitation region, 
its geometry has considerable influence on the antenna impedance or admittance, although 
it practically does not influence the antenna radiation pattern. Therefore precise modeling 
of the antenna excitation region is necessary for obtaining accurate theoretical values of 
the antenna impedance or admittance." 

From the preceding discussion it is clear that rigorous analytical formulation of the 
microstrip antenna problem alone is not sufficient to guarantee accurate input impedance 
results. Careful consideration of the excitation region appears to be just as important in 
input impedance calculations. Moreover, this aspect of the input impedance problem is far 
from trivial, leading certain investigators to conclude that [2] "the most difficult aspect of 
the patch radiator problem is the modeling of the feed. For one thing, many different 
types of feeds are possible: probe feeds, microstrip line feeds, proximity feeds, slot or 
aperture feeds, etc. Some of these feeds are very useful from a practical viewpoint, but 
are very difficult to treat analytically." There is little doubt that the feed modeling prob- 
lem merits a great deal of attention. It will be the primary task of this thesis to study 
and model coaxial probe and microstrip transmission line feeds in the context of the 
overall microstrip antenna analysis. 

The need for theoretical analysis of probe-fed microstrip antennas with electrically 
thick substrates is motivated by several major factors. Among these is the fact that 
microstrip antennas are currently being considered for use in millimeter-wave systems. 
The substrates proposed for such applications often have high relative dielectric constants 
and, hence, appear electrically thick. The necessity for greater band with is another major 
reason for studying thick substrate microstrip antennas. Consequently this problem, and 
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particularly the input impedance aspect of it, has received considerable attention in recent 
years. Several methods, varying in accuracy and computational effort requirements, have 
been proposed and used to calculate the input impedance and other parameters of probe 
fed microstrip antennas. The simpler of these methods is exemplified by the cavity 
model, in which the probe feed is replaced by a uniform current ribbon of equivalent 
width. The cavity model has been successfully used to calculate the input impedance of 
microstrip antennas with substrates ranging in thickness from 0.005 k e to 0.02 X e (\ e - 
wavelength in the dielectric) [3,4]. For thicker substrates, however, the cavity model, as 
well as other simpler models [22], have been abandoned in favor of the more rigorous 
moment method computations of the antenna patch currents. In references [5,6,21] 
moment method formulations were used to solve for currents on the rectangular patch 
fed by a current filament of 1 A. The effect of the probe was then simulated by adding a 
correction term to the calculated input impedance. In reference [7], a circular patch 
antenna fed by a uniform cylindrical current filament of finite radius was investigated. 
After removing the singular part of the current due to the feed, the patch currents were 
computed by using the vector Hankel transform and Galerkin’s technique. Theoretical 
and experimental input impedance plots in [5-9] were compared for substrate thickness (t) 
ranging from 0.01 X € to 0.035 \ 6 and no verification was given for t greater than 0.035 A e . 
All of the formulations outlined above did not consider the actual feed structure in the 
analysis, replacing it instead by some idealized source. 

An alternate approach to the problem of the probe fed circular microstrip antenna is 
presented in Chapter 2 of this thesis. As was already noted, a proper account of the 
metallic nature of the feed probe in the formulation could be very important, if not more, 
in the input impedance calculations. Therefore, instead of solving the open microstrip 
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structure problem using an idealized feed model as was done in [5-8], a solution satisfying 
the boundary conditions on the probe is obtained, utilizing the f ramework of the cavity 
model [3,4] so that the simplicity of using modal expansion in computation can be 
exploited. As in [4] all the losses, including the radiation, surface wave and copper losses 
are lumped into an effective loss tangent 8 eff . An edge extension formula [10] is used to 
account for the fringing of the fields at the edges of the patch. The analysis presented 
applies to the arbitrary location of the probe under the patch except for the uninteresting 
case of the probe very close to the edge. 

In addition to the probe feed, a commonly employed method of exciting the printed 
patch antenna is by means of a microstrip transmission line. There exist several canonical 
feeding arrangements utilizing the microstrip line, namely the direct edge connection, 
capacitive edge feed, proximity feed, etc. Majority of past investigations of patch anten- 
nas fed by the microstrip transmission line were limited to the direct edge connection case 
[2,5, 6, 9]. Moreover, in many instances no distinction was made in the treatment of probe 
and microstrip line feeds, as the ideal current filament was used to model both types of 
feeds [2-4,6]. More realistic approaches to the problem of the proximity fed microstrip 
dipole were presented first by Oltman and Huebner [31], and recently in a more rigorous 
fashion by Katehi and Alexopoulos [27]. 

A method for computing the input impedance of a circular patch antenna excited by 
a microstrip transmission line is presented in Chapter 3. A distinguishing feature of the 
present analysis is the incorporation into the formulation of the boundary conditions on 
the feed. This kind of rigorous treatment has several advantages, namely it yields accu- 
rate input impedance results, and is applicable to a wide variety of feeding arrangements 
utilizing the microstrip transmission line. The problem is formulated in the Fourier 
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Transform domain. Moment method is applied to solve the integral equation for the 
currents on the patch and portion of the feed line. Considerable attention is given to 
acceleration of convergence of the Sommerf eld-type integrals which arise in the course of 
analysis. 

The concluding chapter presents a summary and a discussion of the obtained results, 
pointing out the extent to which the original goals of this investigation were accomplished, 
and leading to suggestions for possible improvements and extensions of the present work. 
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CHAPTER 2 

INPUT IMPEDANCE OF A PROBE FED 
CIRCULAR MICROSTRIP ANTENNA WITH THICK SUBSTRATE 

A method of computing the input impedance for the probe fed circular microstrip 
antenna with thick dielectric substrate is presented in this chapter. Utilizing the frame- 
work of the cavity model, the fields under the microstrip patch are expanded in a set of 
modes satisfying the boundary conditions on the eccentrically located probe, as well as on 
the cavity magnetic wall. A mode-matching technique is used to solve for the electric 
field at the junction between the cavity and the coaxial feed cable. The reflection 
coefficient of the TEM mode incident in the coaxial cable is determined, from which the 
input impedance of the antenna is computed. Measured data are presented to verify the 
theoretical calculations. Results of the computation of various losses as a function of sub- 
strate thickness for the circular printed antenna are also included. 

In Section 2.1, a set of orthonormal modes for the circular waveguide with perfect 
magnetic (PMC) outer conductor and a perfect electric (PEC) inner conductor is con- 
structed. The cutoff wavenumbers for this configuration are defined as the zeros of an 
infinite order determinant, whose elements are transcendental functions of the 
wavenumber. However, a truncated determinant can be used to compute the lowest 
cutoff frequencies with excellent precision [11,12]. The set of transverse modes, 
corresponding to the calculated cutoff wavenumbers, is found next. Each mode satisfies 
the boundary conditions on the outer (PMC) and inner (PEC) conductors. A procedure 
for normalizing the modes over the eccentric annular cross-section of the waveguide is 


outlined. 
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In Section 2.II a mode-matching technique is used to solve for the electric field in the 
annular aperture of the coaxial feed cable. The reflection coefficient of the TEM mode 
incident in the coaxial cable is determined, from which the input impedance of the 
antenna is computed. 

Section 2.III contains a discussion of the results. Theoretical and experimental input 
impedance plots are compared for various values of substrate thicknesses, ranging from t 
- 0.03 \ € to t - 0.1 k e . Finally, various losses for the circular microstrip patch versus its 
substrate thickness are presented. 

2.1. Modes for the Eccentric Annular Waveguide 

According to the cavity model, the circular microstrip antenna can be analyzed by 
assuming that the fields under the patch are essentially those of a magnetic wall cavity. 
The cavity is fed by a coaxial cable. The inner conductor of the cable spans the substrate 
and makes a contact with the patch. Figure 2.1(a) illustrates the resulting configuration. 
A transverse cross-section of the cavity, including the relevant coordinate systems, is 
shown in Fig. 2.1(b). The aim of this section is to construct a set of orthonormal modes 
for the waveguide having the eccentric annular cross-section of Fig. 2.1(b). The task can 
be simplified somewhat by making use of certain inherent symmetries. Since the TEM 
mode incident in the coaxial cable has no angular variation, it is expected that only even 
modes will be excited in the cavity. Moreover, it can be shown that the probe current 
excites predominantly the z component of electric field, and consequently the TM Z modes, 
in the cavity. TE Z modes will not be excited as strongly as the TM Z modes, and will be 
neglected henceforth. 
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2.I.A. Cutoff wavenumbers and modes for the eccentric annular waveguide 

Computation of cutoff frequencies and modes for eccentric annular domains, such as 
that in Fig. 2.1(b), has been the subject of several papers [11-13], on which the first part 
of this section is based. The transverse TM Z modal fields can be derived in the following 
manner [14] 


e q = -- 


v A 


h q = z x e q 


(2.1a) 


( Vt 2 + k 2 ) * q = 0 , =0 on C G , tf q = 0 on Cj (2.1b) 

where a single-integer index set {q} is used to enumerate the modes, and C 0 , Cj denote the 
outer and inner conductor boundaries, respectively. A general solution of Helmholtz’s Eq. 
(2.1b) in the (p, d>) coordinate system of Fig. 2.1(b) can be written as 

tf q (p, </>)= £ C mq ff4 m) Ck q p)cos(rn</») (2.2) 


where {C mq } m=0 is a set of undetermined coefficients and 2? n (m) (k q p) is a combination of 
Bessel functions of first and second kind J n (*) and Y n (*), respectively. Although other 
combinations are possible, the particular one used here is given by 


a„ < ” > (k qP ) = j,(k,p)- 


J m (k,a) 


Y„Ck,p) 


(2.3) 


Note that B^ ml (k q a) — 0. which in torn implies that the boundary condition ip. = 0 on C ( 
is automatically accounted for by the expression in Eq. (2.2). However, the boundary 
condition on C 0 still remains to be satisfied. This is accomplished by using the addition 
theorem for Bessel functions [15], to express tf q in the (p\ <t>') coordinate system as fol- 


lows 
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%(p, <*>') = £ Cjnq £ E(m,n,e) 3 n (m) (k q p')cos(ndO , p > e (2.4) 

m=0 n=0 

where 

E(m.n.e) = |j B _ m (k q e) + (-l) m J n+m (k q e) ] ; y n = \ 

and e is the distance separating the origins of the primed and the unprimed coordinate sys- 
tems. Application of the boundary condition 00 q /Qp' = 0 at p = c and use of orthogonal- 
ity properties for trigonometric functions yield an infinite set of homogeneous equations 

oo 

with {C mq } m=0 as the unknowns. In order for this system of equations to have a non- 
trivial solution, its determinant must vanish. Denoting the determinant by D^Ckq) and 
its elements by d nm (k q ), the equation which the cutoff wavenumbers k q must satisfy is 
given by 

Doo(k q ) = 0 (2.5a) 

with 

d nm (k q ) = E(m,n,e) 5' n (m) (k q c), m,n = 0,1,2... (2.5b) 

where the prime indicates differentiation with respect to the argument. A truncated, Mth 
order restriction of D 0O (k q ), denoted by D M (kq), is used in the actual computation. The 
discussion of numerical techniques used in computing the k q ’s is relegated to Section 2. III. 

M 

To complete the construction of the mode, the set of coefficients {C mq } m=0 must be 
found. From the numerical standpoint, this is a problem of solving a homogeneous sys- 
tem of equations with a vanishing determinant. Since the equations are dependent, 

M 

can be set to equal one, and the resulting linear system solved for {C mq } m=1 . 


11 


2.I.B. Normalization of modes in waveguides with eccentric annular cross-section 


To be useful in subsequent calculations, the modes must be normalized over the 
cross-sectional area (5) of the guide. Mathematically, this is equivalent to finding a con- 
stant N q , such that the following condition is satisfied 



• e q ds = 1 


( 2 . 6 ) 


where e_ is defined as e = _JL. 
4 N q 


Since the modes satisfy homogeneous boundary conditions on C 0 + Cj, Green’s 
theorem can be used to derive the following expression for N q : 


N q 2 = If e q ' e q d S = J J l/f q 2 ds 


(2.7) 


Evaluation of the integral in (2.7) can, of course, be performed by a straightforward sub- 
stitution of expression (2.2) (or (2.4)) into (2.7). However, it is clear that such a substi- 
tution would result in a double summation of terms containing integrals of the following 
type 

/ J JjS4 m) (k q p) 5, U) (k q p)cos (m <t>) cos (/</>) j ds (2.8) 


Because of the complex nature of the eccentric annular domain of integration, a closed- 
form expression for Eq. (2.8) is not possible in general. Therefore, the straightforward 
approach to evaluation of Eq. (2.7) would suffer from computational inefficiency. 

Fortunately, there exists an ingenious procedure for handling this problem. The 
detailed derivation can be found in [16] and will not be repeated here. The final result 


states 
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f J*' 2is =c,L 


w* a* q _ aX 

6(k q 2 ) 9n q QnQ(k q 2 ) 


dc 


(2.9) 


This expression, after further simplification by using the boundary conditions in Eq. 
(2.1b), becomes 




'/'q ds = jf 




(2.10a) 


Close examination of Eq. (2.10a) reveals several important facts. First and foremost is 
the fact that evaluation of the contour integrals becomes trivial, when appropriate expres- 
sions for xlr q are employed. This means that expressions (2.2) and (2.4) should be used on 
contours Cj and C 0 , respectively. Moreover, i f can be shown using expression (2.2) that 
the derivative of i/> q with respect to k 2 , evaluated at k 2 = k q 2 and p = a, is identically zero. 
This implies that the contour integral on C t must also vanish. 


As noted in the previous paragraph, evaluation of the remaining contour integral on 
Cq is trivial. The result is in the form of a triple summation, which can be written as fol- 
lows 


2 °° °° 

N 2 = 770 


I L c* c„ z A 


* i=0 ]=0 k=0 z 


I 

(e/c) Jj'n.jtkqe) + (— l)‘J' k+i (k q e) 


3' k (i) (k q c) + 


J k -i(k q e) + (— l)'J k+i (k q e) ] 


B"«K k q c) + _L_ ntkqC) 
q TT k q c Y 2 ( k q a) 


x Jj k _j(k q e) + (— l) j J k+ j(k q e) |s k (j) (k qC ) 


(2.10b) 


where 
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S'^>(k,c) = J' k (k,c) - Y' k (k q c) 

ri<K k,c) = r t (k,c) - Y-;(k,c) 

and primes indicate differentiations with respect to argument. Numerical evaluation of 
the summations in Eq. (2.10b) is several orders of magnitude faster than the straightfor- 
ward evaluation of Eq. (2.7) original proposed. 

To recapitulate, a general derivation of modes for an eccentric annular waveguide 
was presented in this Section. A dispersion relation for the cutoff frequencies of such a 
waveguide was given, along with a method for finding the corresponding modes. Finally, 
a very efficient way of normalizing the modes over the cross-sectional domain of the guide 
was presented. 

In the special case wherein the ratio of inner to outer conductor dimensions is small, 
perturbation techniques can be applied to derive analytical expressions for the cutoff 
wavenumbers and the corresponding modes. The procedure is illustrated in Appendix A, 
where annular cross-section waveguides with PEC and PMC outer boundaries of both cir- 
cular and rectangular shape are considered. 

2.II. Computation of Input Impedance 

Consider the magnetic wall cavity fed by a coaxial cable, as illustrated in Fig. 2.1(a). 
It is assumed that a TEM mode is incident upon the annular aperture (A) formed by the 
junction between the coaxial cable and the cavity. In this section an integral equation for 
the tangential electric field E a in the aperture, A : a ^ p < b, is derived. The integral 
equation is solved numerically by the mode-matching technique for the reflection 
coefficient (R) of the TEM mode, from which the input impedance for the circular 
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microstrip antenna is determined. 

2.II.A. Derivation of the integral equation for the fields in the junction aperture A 
The fields in the cavity region with permittivity € t can be represented in terms of the 

oo 

set of orthonormal modes {ei} i=1 derived in Section 2.1. At z - 0 + , the components of the 
field transverse to z can be expressed as 

E t (p,z=0 + ) + £ V^Cp) (2.11a) 

i=i 

-z x H t (p,z=0 + ) = £ V^p) (2.11b) 

i=l 

where Y it i-1,2,... are the modal admittances, at z = 0 + given by 

Y ; = — j cot(kj.jt) ; k 2i = V<<>Vo € i~ k i 2 (2.12) 

^■zi 

Using the orthonormality property of the efs and Eq. (2.11a), the modal coefficients, V j( 
can be expressed as follows 

= E t (p,z=0 + ) • e;(p)ds = ff E t (p,z=0 + ) • ej(p)ds ; i=l,2,... (2.13) 

where the domain of integration is reduced f rom S to A on account of the boundary con- 
dition which requires that E t (p,z=0 + ) vanish on S—A . 

The coaxial waveguide modes are well-known [17], and are briefly presented in Table 
2.1 for future reference. At z=0“, the transverse fields in the coaxial cable can be 

OO 

represented in a complete set of orthonormal modes as follows 

E t (p,z=0") = (1+R )e j(p) + £ Vje'j(p) 


(2.14a) 



MODES FOR THE COAXIAL WAVEGUIDE 



Potential Function 



TEM Xoo = -i2W- 

Moo 


®oo — “V tXoo 


™.' v “ ! = ^^flcosCm*) , m=0,l,2. 


e pm ~ 


^7 t X 0 


TE z even x ttn = - ”^ Qnp) sin(n <*>) , n=l,2,3... 

M 2an 


S' _ ^tXan v ~ 

e an _ A z 

*om 


FaS-K-vmP^ ^aS- K vmP^ m^vm^ 3aS- K viifl^ m {K vm p) 


Gn(*anP) ^n^ K an P"N ^ n^*an a ) Y,(K m p) 


Mqo = \2ir /n(b/a) 


M. = Of,. ^.'1-1/2^ -1 

ivl li>m ^Ym 7rj K vm rrrf TT 1 

Jm(K wm b) 


M = 9— -1/2^ -1 Jn( K an a ) , n 2 

*“ ^ TW7 ‘-CCb? 


Ulna) 2 
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— z x H t (p,z=CT) = (l-R)yie j(p) - £ v^-e fp) 

j=2 


where y^’s are the modal admittances given by 
yj = V * 2/^0 • TEM mode (j=l) 


(2.14b) 




Z] 


, TM_ modes 


2) 

“Mo 


, TE- modes 


(2.15) 


and 


K zj = -v/“ 2 Mo« 2 “^j 2 * 

The e 2 is the medium permittivity in the coaxial waveguide and the k } ’s are the cutoff 
wavenumber of the coaxial waveguide, enumerated by the single-integer index set {j}. 
Again, since the infinite basis {e j} is orthonormal, it is possible to define modal coefficients 
Vj, j-1,2,... as 

v j = // E t (p,z=0-) • e j(p)ds (2.16) 

where Vj - 1+R. 

Following the theory outlined in [19], expressions (2.11b) and (2.14b) are equated in 
order to satisfy the continuity conditions of the fields on A . Substitution of Eqs. (2.13) 
and (2.16) into the resulting equation and interchange of order of summation and integra- 
tion [19] yield 


2y,e ,(p) = // 


£ yje j(p')e j(p) + £ Y&Cp'KjCp) 


U =1 


i=i 


Ejip'W 


(2.17) 


Equation (2.17) is the sought integral equation for the electric field E a in the aperture A . 
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2.II.B. Solution of the integral equation by the mode-matching technique 

The method for reducing the integral equation in (2.17) to a set of linear algebraic 

equations is well-known [19]. To briefly summarize, the held E a (p') is chosen to be 

K 

approximated by a finite set of modes |e n } n=1 , i.e., 

E a (p')~ £ v„e n (p') (2.18) 

n=l 

and the so-called intermodal coupling coefficient is defined in the following fashion: 

r ij = // * e j cis (2.19) 

A 

Computation of I» is in many instances facilitated by the use of Green’s theorems and 
other vector calculus identities, in the manner demonstrated in Appendix B. The final 
expressions for the inner products between the coaxial cable modes, summarized in Table 
2.1, and the cavity modes defined in Section 2.1 are tabulated below. 

Following [19], expression (2.18) is substituted into (2.17) and moments are taken 
with e m (p), m«l,2,...,K to obtain a linear system of equations given by 

IGIrxkMrxi = WkxI 

G m „ = y m 8 mn +£ YiI im I in (2.20) 

i=i 


2yi8 lm 


8 m „ = 


1 , m=n 
10 , m^n 


The KxR matrix [G] can be inverted to determine the modal coefficients Vj, j«l,...,K. 
The reflection coefficient R and the normalized input impedance Z in are then determined by 
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Table 2.II 


1 II 

| INNER PRODUCTS I i} = /JVe'jds jj 

1 * II 

i ii 

Mode types 
cavity 

coax 

— — — — 1 - - ■■■ ■ ■ 1 n 

TM Z ven — e q 

wM 

N q M»k q C ^o W >(k q b) 

TM Z ven - e m 

1 

T —4 1 (kqlK. vw ) n n (m)ft Kt 

M**) ym NqM]>m U4 _ k 2) jJ^bT C “> ” 1 ’ 

1 

! TE z even — e an 

Iw ‘“ ) N q M 2 „ kJ^Fb c ”> fl " (k ’ b) 

1 

1 

1 

1 

Note: Definitions of all the quantities appearing in this Table can be 
found in Table 2.1, Eqs. (2.2,2.3,2.10b), Figs. 2.1(a-b). 


R = Vj — 1 , and Z in = - (2.21) 

The coefficients Vj can be calculated by using the following relation 
K 

Vi = £ IjjVj i=l,2,...,L (2.22) 

j=i 

Once the V{s are available, the total electric field in the cavity can be written as [14] 
E c (p,z) = £ Vj (z)e;(p) + - r J_ £ IjCz) V t ' CjCp) (2.23) 

i=i i=1 


where 
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V,(z) = V, 


sin[k zi (t— z)] 
sin(k 2j tj 


Ij(z) = j 


dVj 
TT d z 


2.II.C. Discussion of the "effective dielectric loss" tangent calculation 

Thus far, the microstrip antenna has been modeled by an ideal closed cavity. How- 
ever, the antenna proper exhibits various loss mechanisms, including radiation and 
surf ace- wave excitation, and a method for incorporating these into the model must be 
found. One possible method is to lump the radiation, surf ace- wave, copper and dielectric 
losses into an "effective dielectric loss" tangent and replace in Eq. (2.12) by 
1— j 8 eff ), as suggested in [4,18,20], The quantity 8 eff is defined by the following 
expression 

Seff = 2 o) t W7 (2.24) 


where P - P rad + P sw + P cu + P d is the total power lost due to radiation, surf ace- wave 
excitation, imperfectly conducting metallic portions of the cavity and in the dielectric. 
W E is the time-averaged stored electric energy at the resonant frequency to r . 


2.II.D. Evaluation of radiated and surface-wave power 

Several methods exist for computing the radiated and surface wave power. The 
method used here assumes that a magnetic current M = E c xn | ^ is flowing on the outer 
PMC wall of the cavity, where n is the unit normal [3] outward from the cavity wall. 
The fields produced by M are found using the Fourier integral representation of the mag- 
netic current Green’s function for the grounded slab [14]. The surface-wave fields are 
found by evaluating the pole contributions of the Fourier representation integrals. The 
method of steepest descent is used to determine the space wave contribution of M. 
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As an example, consider the Green function for a transverse magnetic current source 
embedded in a PEC-backed dielectric slab of thickness t, and with constitutive parameters 
= ^s^i, yti 0 . The space- and surf ace- wave fields excited by the time-harmonic point 
source M = f 8(p— p')8(z— z), r being an arbitrary unit vector transverse to the z axis, and 
(p',z') denoting the source location, are given by the following set of expressions: 

Radiated or space-wave fields 


E rad (r,r') = e^jk,, x (jk 0 x z) n' rad (r,r') - j <op 0 (jk 0 x z)II" rad (r,r') 
H rad (r,r') = j<o€j (jk 0 x z) II' rad (r,r') + jk 0 x (jk 0 x z) II" rad Cr,r') 


(2.25a) 


j®«i n' rad (r,r') = f x z • V't G' rad (r,r') 
jo>p 0 n" rad (r,r') = -J_ f Jjt • V't G" iad (r,r') 


(2.25b) 


C^Cr.r'J = |i"£« Y' rwi (z^')g 1 (k 0 rV it ' p, * inflc0 * < * - *' ) 

G" tia (r.r-) = 8,0V) ?*****-< 


) 


(2.25c) 


where 


Y- (7 = o>€ 1 c os(k zl z' )e jklot 

rad ’ lk zo € rl cos(k zl t) + jk zl sin(k zl t)J 

7 - f m _ jo>/u 0 sin(k zl z')e jkwt 

radU,Z J Lk zl cbs(k zl t) + jk zl sin(k zl t)J 

r - . -jV 

k zo = k o cos0, k zl = k 0 V e^-sin 2 #, g^r) = — 

477T 
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V't is the transverse "del" operator, applied to the source coordinates. 


Surface-wave fields (TM Z modes ) 


E sw (r,r') 


*1 

dz) 


V x V x z IT sw (r,r') 


H sw (r,r') = V x z II' sw (r,r') 


dz) = 


«*,. t<z 
€j, 0<Z<t 


(2.26a) 


j“«i n' sw (r,r') = f x z • V't G' sw (r,r') 


(2.26b) 


G' (r r') = 1 Y' 

^ sw vi ** > 25T 1 5 


,(z,z')g 2 (k c p) e 


jl^p'cosC^-^') 


(2.26c) 


where 


Y' sw (z,z') 


WPJ 


cos(/3 1 z')cos(/3 1 z) 

cosC^t) 

cos^z'le t} 


. , 0<z^t 
, t<z 


D E (3 0 ,3x) = sin(3 x t) 


k 0 2 (^i-D ^ tffi + ftfej) 

3 0 2 3i Ml«rl 


k 0 = Vk c 2 -k 0 2 , /8, = VkJ^k?; g 2 (k c p) = J2hrt~p 


dispersion relation for k c : (/3jt) = (/3 0 t) ^xCOtl /3x t) 
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Surface -wave fields ( TE Z modes ) 


E sw (r,r') = — ju>fx 0 V x z n" sw (r,r') 

H sw (r,r') = V x V x n" sw (r,r') 

joj/z 0 n" sw (r,r') = _L_ r ^ • V' t G" sw (r,r') 

0 2 

G" sw (r.r’) = i Z\Jz, z) g,(k c p) 
where 


(2.27a) 


(2.27b) 


(2.27c) 


Z" (z,z) 


W/OT 


sin(j3iz’ ^inC^z) 
sin(/3 1 t) 

sin(j8z')e ^ 0<z l) 


, 0<z<t 
, t <z 


D h ( 3 0 ,/3 1 ) = — cos(3it) 


(l+^tX^+fl, 2 ) 

§73; 


dispersion relation for k c : (3 0 t) = — (3it)cot(0it) 

The power components P rad and P sw are calculated by integrating the corresponding 
Poynting vectors over the appropriate regions of space. Results of the numerical evalua- 
tion of various components of P for a specific antenna are presented in the following sec- 


tion. 
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2.III. Numerical and Experimental Results 

Numerical results based upon previously performed theoretical analyses are 
presented in this section. Theoretical input impedance data are verified experimentally in 
all cases. 

In Section 2.1 a dispersion relation for the cutoff frequencies of a circular PMC 

waveguide with an eccentric inner PEC conductor was derived. The dispersion relation 

was given in the form of an infinitely dimensional characteristic determinant D^Ckq). 

The cutoff wavenumbers were defined as the zeros of this determinant. In actual compu- 

12 

tations, a truncated determinant of order M was used. Calculations of {k q } q=1 were per- 
formed for various ratios of a/c. It was found that the maximum discrepancy between 

5 

the results for M-10 and M=»15 is confined to within 1%. A graph of <k q c} q=1 is shown in 
Fig. 2.2. The ratio of inner to outer conductor radii for this particular case is a/c - 0.017. 
Note that as e goes to zero, the cutoff wavenumbers approach the well-known values of 
k q c for the concentric annular waveguide. In fact, these values can be utilized as initial 
estimates for the numerical root searching routines. 

In Section 2.II a method of computing the input impedance Z in for the circular 
microstrip antenna was presented. A mode-matching technique was used to solve for the 
reflection coefficient R of the TEM mode in the coaxial cable, f rom which Z in was calcu- 
lated. In all cases considered, the cavity radius c was much larger than the outer radius b 
of the coaxial cable. Consequently, the number of modes required to approximate E a in 
the coaxial cable region was found to be considerably less than the corresponding number 
in the cavity region. Approximately eight to ten cavity modes were required to obtain 
numerical convergence. In contrast, only one or two coaxial cable modes were found to 
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give very accurate results, avoiding the need for inversion of large matrices. In fact, the 
problem could have been formulated from the outset by assuming a TEM field distribu- 
tion in the aperture A and applying the variational method [18] to determine the input 
impedance. 

In Figs. 2.3-2.5, Smith Chart plots of theoretically computed values of input 
impedance are compared with measurements for circular microstrip antennas of various 
sizes. The substrate material used in all cases was Rexolite 2200 with relative permit- 
tivity of 2.62 and a loss tangent of approximately 0.001. The input impedance loci were 
computed for two feed positions in each case. The curves numbered 1 and 2 in Figs. 2.3- 
2.5 correspond to feed locations of 7.5 mm and 10.5 mm, respectively, from the center of 
the patch. As noted previously, an edge extension formula [10] was used to account for 
the fringing field effect. As a result, the shift between the theoretical and experimental 
resonant frequencies is eliminated. 

The results presented in Figs. 2.3-2.5 are for antennas with substrates ranging in 
thickness from 1.6 mm to 4.7 mm. In terms of the dielectric wavelength \ € , this is 
approximately equivalent to having t range f rom 0.03 \ € to 0. 1 \ e . As the substrate thick- 
ness increases, analysis of input impedance data reveals several interesting trends. The 
expected shift of the impedance locus to the inductive side of the chart, as t increases, is 
clearly observed in Figs. 2.3-2.5. The shift is due to the fact that for large t the self- 
inductance of the probe constitutes a significant part of the overall antenna impedance. It 
is also observed that the bandwidth of the antenna markedly increases with substrate 
thickness. The band widths for SWR < 3 and the Q factors for the three cases shown in 
Figs. 2.3-2.5 are given in Table 2.III, along with the relevant antenna dimensions. 
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X X COMPUTED LOCUS 

INCREMENT: O.iGHz 

(INCREASING FREQUENCY IS CLOCKWISE) 

Figure 2.3. Measured and computed impedance loci: c - 14.1 mm, t - 1.6 mm, ej - 7.5 
mm, e 2 - 10.5 mm, - 2.62. 
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x x COMPUTED LOCUS 

INCREMENT: O.l GHz 

(INCREASING FREQUENCY IS CLOCKWISE) 

Figure 2.4. Measured and computed impedance loci: c - 13.5 mm, t - 3.2 mm, ej - 7.6 
mm, e 2 - 10.5 mm, - 2.62. 
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MEASURED LOCUS 


x x COMPUTED LOCUS 

INCREMENT: O.l GHz 

(INCREASING FREQUENCY IS CLOCKWISE) 

Figure 2.5. Measured and computed impedance loci: c - 13.0 mm, t - 4.7 mm, e t - 7.5 
mm, e 2 - 10.5 mm, = 2.62. 
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Table 2.III 

Estimated Bandwidths and Q Factors 


Antenna dimensions 

0 factor 

Bandwidth (VSWR < 3:1) 

c - 14.1 mm 
t - 1.6 mm 

33 

5% 

c - 13.5 mm 
t - 3.18 mm 

15 

10% 

c - 13.0 mm 
t = 4.8 mm 

9 

18% 


Finally, Fig. 2.6 depicts proportions of various power losses for the circular printed 
antenna of radius 0.25 \ e as a function of substrate thickness. It is seen that while the 
dielectric and copper losses become relatively insignificant the surface wave power 
increases with thickness and as much as 30% of the power can be lost through surface- 
wave excitation. As a result, the radiation efficiency of the antenna can fall below an 
acceptable level. Thus, it appears that for broadside radiation, very thick microstrip 
antennas, although having broad-band impedance characteristics, may not be desirable 
f rom the radiation efficiency point of view. 
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APPENDIX A 


CUTOFF WAVENUMBERS AND MODES FOR ANNULAR 
CROSS-SECTION WAVEGUIDE WITH ECCENTRIC 
INNER CONDUCTOR OF SMALL RADIUS 

A.I. Introduction 

Computation of cutoff wavenumbers for uniform waveguides with eccentric annular 
cross-section has been the subject of numerous investigations [il-13,35]. 1 Various tech- 
niques suggested for that purpose are for the most part rigorous in nature, and require 
considerable numerical analysis. A relatively important limiting case which has not 
received much attention, and yet is of practical interest, is the annular waveguide with 
small ratio of inner to outer conductor dimensions. Configurations of this type arise, for 
example in the analysis of cavities excited by thin probes [36]. 

The purpose of this appendix is to derive approximate analytical expressions for the 
cutoff wavenumbers and the corresponding modal wavef unctions for annular waveguides 
with small ratios of inner to outer conductor dimensions. Derivation of the results is 
based to a large extent on the perturbation technique outlined in the investigation entitled 
"Distortion of Standing Wave by a Strip," found in reference [16]. 

Circular waveguide with eccentrically located inner conductor of small radius is 
analyzed in Section A.II. Cutoff wavenumbers computed therein are compared with data 
numerically generated from a more rigorous formulation 111]. 

Section A. Ill considers the case of rectangular waveguide with arbitrarily located 
small inner conductor. Again, the cutoff wavenumbers and the corresponding modal 
wavef unctions are derived. 

1 Only a partial list of relevant publications is referenced. 
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A.II. Circular Waveguide with an Eccentric Inner Conductor of all Radius 

Approximate formulas for the cutoff wavenumbers and the corresponding modes in a 
circular waveguide with eccentric inner conductor of small radius are derived in this sec- 
tion. A cross-section of the waveguide, with the relevant dimensions indicated, is shown 
in Figure Al. Analysis of the TM modes is presented in detail. The results for the TE 
case can be obtained in an analogous manner. 

The TM modes in a uniform waveguide can be derived from a scalar potential func- 
tion Ms vm satisfying the wave equation 

( V t 2 + kjmM/m = 0 . * (Al) 

and subject to appropriate boundary conditions [14]. In the absence of the inner conductor 
in a circular waveguide having a PEC outer boundary of radius c, the solution of Eq. (Al) 
is known to be 

= Jm( k pmP')cos(rn0') ; J m (k„ m c) = 0 (A2) 

where J m denotes a Bessel f unction of the first kind of order m. It is intuitively clear that 
the presence of a small inner conductor slightly perturbs the above solution. Derivation 
of the modified solution, denoted by % m , is facilitated by the introduction of several 
hypotheses, detailed discussion of which can be found in [16]. A brief restatement thereof 
follows. 

Initially it is assumed that the standing wave which is expressible as a combina- 
tion of diverging and converging cylindrical waves, is incident upon the inner conductor 
and is scattered by it. The effect of the scatterer on the incident waves is computed as 
though these waves were out in the unbounded medium, i.e. with the outer wall boun- 
dary C Q at infinity. It is hypothesized that the value of ? vin at the inner conductor 
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boundary C t should not be very sensitive to the presence or absence of the outer wall, so 
long as the inner conductor is small and is not situated close to C 0 . Consequently, the 
behavior of on Cj is known fairly accurately from the analysis performed with C Q at 
infinity. The behavior of the modified wavef unction in the rest of the waveguide 
cross-section can be found with the aid of the following equation 

^mCpo’^o) - -/ (Ps A ) GCp 0 ,d> 0 |p s A )dc , (A3) 

c, 6 n 

derived using Green’s theorem. The subscripts "o" and "s" in the above equation refer to 

t ». / » 

the observation and source coordinates, respectively. G(p 0 ,</> o |p s ,<2> s ) denotes the scalar 
Green function satisfying the appropriate boundary condition on C Q . In the case of a PEC 
outer boundary, application of standard analytical techniques yields the following expres- 
sion for G(p 0 ,d» 0 |p s ,<£, ) 


G(p o '.0 o |p s ,d> s ) = Z Z 

n=0 f=l 


^qnPo nPs ) 

N/ n (k 1 2 -k/J 


cos n{<t> Q — <b . ) 


where 

1/2 

Jn+AfnC) CA4b) 

1 , m=0 

y n = 

2 , m^O 

In accordance with the procedure outlined in preceding paragraphs, let it be assumed 
that the combination of cylindrical waves, given by xlt vm , is incident upon the inner con- 
ductor. The conductor is centered at the origin of the coordinate system (p,<t>), shown in 
Fig. Al. The addition theorem for Bessel functions [15] is used to transform (Eq. 
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(A2)) from the (p',d>') into the (p,<t>) coordinate system. The resulting expression is given 

by 


\lr vm (p,<t>) = £ E(m,p,e)J p (k um p) cos(pd>) (A5a) 

p=o 


where 


E(m,p,e) = |j m _ p (k wm e) + (-D^+ptk^e) ] 

The scattered field can be expressed as a sum of outgoing cylindrical waves 2 


(A5b) 


^ = £ A p Hp 2) ( k wm p )cos( pd>) 

p— 0 


CA6) 


where Hp r) is the Hankel function of the second kind of order p. The unknown 
coefficients A p are found by enforcing the boundary condition requiring + rl? to 

vanish on Cj(p=a). The incident-plus-scattered wave solution and its normal derivative 
evaluated on Q are given by 


= £ E(m,p,e) 

p=0 


ll fir ^ - J P ( - k «'«n a) .u( 2)rt. 

lUKmP) H^'k^aJ p vmP 


cos(pd>) 


(A7) 


ay 


a^ vm (p,d>) . 


an 


’ t - t 1 '- - - Hr £ E( - e) Ti-cbr • 


jra 


CA8) 


where the Wronskian of Bessel’s equation was used in deriving Eq. (A8). In the limit, as 
the radius a goes to zero, Eq. (A8) can be approximated by the following expression 


a^vm i 2 E(m,0,e) __ 2 
Qn |C| Ira fCk vm aj 7ra f(.k„ m aj 

where f(k„ m a) is the small argument expansion of jH 0 (2) , given by [15] 


e- 1 " 1 time dependence is assumed throughout the thesis. 


(A9a) 


2. 
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nu = - 


log 


Km* 


+ y 


(A9b) 


and y - 0.5772 is Euler’s constant. 


At this point, Eqs. (A4.A9) can be substituted into Eq. (A3) and the prescribed 
integration performed. The most expedient way to perform the integration is by 
transforming all the functions involving the variables of integration into the (p,<*>) coordi- 
nate system. In the case of the Green function G(p 0 ,d> 0 |p s ,<5 S ), the required transforma- 
tion is facilitated by the use of the addition theorem for Bessel functions. The transfor- 
mation is applied to source coordinate (p s ,<5 S ) and yields the following expression 


G(p o ,d> 0 |p s A )|c, = Gtp'.dfJp.dOl^ = £ £ 

n=0 f=l 


Jn(kfnP') 

N f i,u 2 -kjy 


• £ E(n,q,e) J q (k fn a) [cos(nd>') cos(qd>) + sin(nd>' )sin(qd>)] (A10) 

q=0 

where the subscript "Of* is unambiguously omitted from p',d>. Equation (A3) can now be 
rewritten as 


Itt 

v um {PA)~- / 




9p 


G(p',d>'|p,$) 


|_ a • (add>). 


(All) 


straightforward evaluation of which yields 

\ir (n eh'} — 4J m (k„ m e) * ” J n (k| n e)J n (k^ n p ) , ■> 

WP 3 = “ - f ( k ’ I I M 17 ' ul , h ’ Jo^{ n a)cos(n <t> ) (A 12) 

HK vm aj „=o f=1 N^fk^— k£,) 

As it stands, Eq. (A 12) is not correct unless the parameter k is replaced by the 
modified cutoff wavenumber K vm . Noting the fact that xlJ vm and satisfy the wave 
equation, with wavenumbers k„ m and K„ m , respectively, and using Green’s theorem, it can 
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be shown that K pm satisfies the following formula [16] 


K 2 — i, 2 _ 

“■urn 



(A13) 


where S A is the annular region bounded by C 0 and Cj. Using the hypothesis that does 
not significantly differ from \lt vm over most of S A , the integral in the denominator of Eq. 
(A 13) can be approximated in the following manner 

*vm ds 35 fjl d S = N v 2 m ( A14) 


where S c is the circular area bounded by C Q , and N vm is given by Eq. (A4b). Substitution 
of Eqs. (A5a) and (A9) into Eq. (A13) yields the following result 


K 


2 


4 Jm(k„ m e) J 0 (^ m a) 
N y 2 m f(k vm a) 


(A15) 


Note, the first term in the above expression for K„ 2 m is just the square of the original cutoff 
wavenumber, whereas the second term constitutes a small perturbation thereof. 

To complete the derivation of a substitution of K vm for k in Eq. (A12) is made, 
and the result simplified using the fact that 


V 2 _ u 2 

rv vm 1 '-fn 


k 2 — k r" 

'n'm 


, |n 5* I'm 


4J"(K, m e)J 0 (k 

vm = vm 

N„ 2 m f(k„ m a) 


(A16) 


the final expression for can be written as follows 






4Jmtk,„e) ^ 
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OO OO 

X I I 

n=0 f=l 
fn * I'm 


Jn(k^ n e)J n (k (n p) 

N/„ (k, 2 m - k/ n ) 


J 0 (k| n a) cos(n <£') 


(A17) 


where in addition to the original wavefunction xlt vm there is a small term representing the 
perturbation by the inner conductor. 

To verify the accuracy of the presented solution, a comparison is made between the 
results computed from Eq. (A 15) and those obtained by the more rigorous technique of 
reference [ill. Lowest eight cutoff wavenumbers are plotted in Figs. A2, A3, as functions 
of the eccentricity e/c, for a/c ratios .01 and .03, respectively. The solid curves represent 
the data computed from Eq. (A15) and the markers are placed at points computed by the 
alternate method [11]. The agreement between the separately computed sets of data 
appears to be very good for the cases considered. 

Clearly, the pertubration procedure outlined in this section is equaly applicable to 
the case where the outer boundary C 0 is a Perfect Magnetic Conductor (PMC). In fact the 
derivation of K„ m and in that case is identical to the one presented above, and the 
results are formally given by Eqs. (A15) and (A17) with N„ m replaced by 

1/2 

(A18) 

Note, the results do not apply for the lowest cutoff wavenumber, which in the 
absence of the inner conductor is equal to zero. However, this does not present a major 
problem, since it can be readily computed by other methods. 

Fig. A4 plots the modified cutoff wavenumbers for a waveguide with a PMC circular 
outer conductor and a PEC inner conductor. The ratio of the radii is a/c - .022. The solid 
curves are computed from Eq. (A 15) with N vm given by Eq. (A 18). The markers are 


N„ m = 


yc* 

Ym 


1 - 


m 


(k„ m c ) 2 


Jm(k vm c) 
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placed at data points obtained by the method outlined in Section 2.1. 


A.III. Rectangular Waveguide with a Cylindrical Inner Conductor of Small Radius 

Analysis of the circular waveguide with an eccentrically located inner conduction of 
small radius was presented in the preceding section. It should be noted, however, that the 
general technique used in Section A.II to derive the cutoff wavenumbers and the 
corresponding modes is applicable to other waveguide cross-section geometries. In this 
section it will be used to solve for the TM modes in a rectangular waveguide with inner 
conductor of small radius. Figure A5 illustrates the geometry, along with the relevant 
coordinate systems and dimensions. To avoid repetition and keep the presentation as brief 
as possible, all the assumptions and hypothesis stated in the preceding section will be 
omitted here. Only the most essential steps in the derivation will be shown, since it is 
completely analogous to the one found in Section A.II. 

In the absence of the inner conductor, the "mn"th TM mode in a waveguide with a 
rectangular PEC outer boundary can be determined from the scalar potential function: 

mn = sin(k xm x')sin(k yn y'); (A19) 


where 


mf 

k = 

njr 

c 

* **yn 

b 


The wavefunction \U mn is expressible as a combination of four plane waves. It is assumed 
that these waves are incident upon the inner conductor. The scattering problem of plane 
wave by a circular cylinder is well documented [16] and its solution need not be recon- 
sidered here. The total incident-plus-scattered wave 'P mn = 'Jt mn + tf*. satisfying the 
homogeneous Dirichiet boundary condition on Cj, is found by superposition of the 
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solutions for the four plane waves constituting \lr mD . As in Section A.II, the quantity of 


interest is the value of 


d* 


mn 


9 n 


k- 


Straightforward analysis shows that for small values of 


d^mn , ^ 2 . sin(k xm x c )sin(k yn y c ) 
0n * C| ira f(k„ m a) 


(A20) 


where k mn = Vk xm +k yn and a , x c , and y c are as shown in Fig. 5. 

The form of the modified wavef unction ^ mn at arbitrary points within the guide 
cross-section is found with the aid of Eq. (A3). For the case of rectangular waveguide 
being considered here, the Green function appearing in Eq. (A3) is given by 


G(x 0 .y 0 |x s ,y s )= £ £ 
p=i q=i 


sin(k xp x 0 )sin(k yq y 0 )sin(k xp x s )sin(k yq y s ) 


N 2 (k 2 — k£) 


(A21a) 


where 



(A21b) 


Integration prescribed in Eq. (A3) is facilitated by several transformations of Eq. 


t * 

(A21). The source coordinates transformation Xj = x c + pcos<£, y s = y c + psin</> is imple- 

» t 

mented first. This allows the term sin(k xp x s )sin(k yq y s ) in (A21a) to be rewritten as fol- 
lows 


• I 

sin(k xp x s )sin(k yp y s ) = [sin(k xp x c )cos(k xp pcos0) + cos( k xp x c )sin( k xp pcos</> )] x 


x [sin(k yq y c )cos(k yq psin0) + cos(k yq y c )sin(k yq psin0)] (A22) 

Equation (A22) is further expanded in a series of terms involving Bessel functions 
J„(k xp p), J„(k yq p) and angular harmonics cos( sin( v<t>) [15]. The first and only term of 
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the expansion required for subsequent calculations is given by 

sin(k xp x s )sin(k yq y s ) 555 sin(k xp x c )sin(ky q y c )J 0 (k xp p)J 0 (k yq p) (A23) 

Substitution of Eqs. (A20) and (A21a) into Eq. (A3) and subsequent use of Eq. 
(A23), yields the formal expression for the modified wavef unction ¥ mn : 


, T , r > ~ 4 sin(k xm x c )sin(k y c ) » “ sin(k xp x')sin(k y y')sin(k xp x c ) 

^ mn (x ,y ) ~ t • T Y • 


f(k mn a 


I I 

p=l q=l 


N 2 (k 2 -kj) 


sin(k yq y c ) • J 0 (k xp a)J 0 (k yq a) 


(A24) 


As previously noted in Section A.II, Eq. (A24) is not correct until and unless k is 
replaced by the modified wavenumber K mn . Equation (A13) is used to calculate K mn> 
which after some analysis can be shown to satisfy the following formula 


, , 4 sin 2 (k xm x c )sin 2 (k vn yJ , . , 

* k^ n N T f(k a) ' ■ — * J 0 (k xm a)J 0 (k yn a) 


(A25) 


The final form of the modified wavefunction ¥ mn is obtained when Eq. (A25) is sub- 
stituted into Eq. (A24) and the resulting expression simplified in the manner previously 
described in Eq. (A 16): 


*_(*✓) = - 4sinU ^ (k ^ 


f(k nn aJ 


z z 

p=i q=i 

pq * mn 


« sin( k xp x' )sin( k y q y ' )sin( k xp x c )sin( k yq y c ) 

L rrrr. 7 — m Jo^xp 3 ' JotK yq aj 


N^k^-kj 2 ) 


yq a '(A26) 


Figures A6(a,b) plot the variation of K mn with x c , for a rectangular waveguide with 
ratio b/c - 0.5. The results are plotted over the domain 0.5 ^ x c /c ^ .97, y,/c = 0.1, 0.2. 





